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Abstract 

A question is given on the form n(fip — /x Q ) for the volume term of work of 
formation of critical nucleus. Here, n is the number of molecule undergone 
the phase transition, /i denotes the chemical potential, a and (3 represent 
the parent and nucleating phases, respectively. In this paper we concentrate 
phase transition without volume change. We have calculated the volume 
1 _Pj! term in terms of the chemical potential difference /i re — fx eq for this case. 

Here, \x re is the chemical potential of the reservoir and [i eq that at the phase 
transition. We have 

gj: w vol = -^^^ re -^ q yv p 

o ' 

with k denoting the isothermal compressibility, v eq being the molecular vol- 
£D • ume at the phase transition, Vg the volume of the nucleus. 
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1. Introduction 

To calculate the reversible work of formation of the critical nucleus is one 
of the main purposes of the theory of nucleation, because one can predict 
the steady-state nucleation rate J s = J exp(— W* /UbT) through the work of 
formation of the critical nucleus W*, where /c#T is the temperature multiplied 
by Boltzmann's constant. Here, W* = W(R*) is the height of the work of 
formation of critical nucleus with R* being the radius of the critical nucleus. 
We often encounter the following formula, Eq. (p^), or equivalent one: 

W = n(/ip - n a ) + 7A, (1) 

with 7 begin the interfacial tension, A = AnR 2 the area of the interface 
(rigorously speaking, R is the radius of the surface of tension) in textbooks 
such as Refs. |l|, |2|, y|, ]% la, la] as well as research papers such as Refs. p], |g, |9|, 



10J, llll, |12|, |13| . In this formula, one regards \ip — \i a as the chemical potential 
difference between the parent phase (the a phase) and the nucleating phase 
(the (3 phase). One can understand n as the numbers of molecules undergone 
the phase transition from the a phase to (3 phase. It seems to imply that no 
volume change is assumed to be associated with the a-/3 phase transition. 
That is, A(n/i) reduces to nA/i in the case that n is common to both a 
and (3 phases; if n a = V$jv a and np = Vp/vp with Vp = 47ri? 3 /3, one has 
n a = n p{= n) for the case that the molecular volumes v a and vp equals with 
each other. In crystal growth from the melt, the difference between densities 
of the crystal and melt phases is often neglected. Indeed, Eq. (JTJ has been 
used for crystal nucleations in the melt. However, while someone describes 
Eq. ([1]) or equivalent one in phase transitions with small volume change such 
as melt-crystal cases, someone does in vapor-liquid cases. In some literatures, 
implication is unrevealed even by reading between lines. 
The exact form for W given by Gibbs [14| is 

W = -{pp-p a )Vp + 1 A, (2) 

where pp and p a represent the pressures of the respective phases. Rigorously 
speaking, pp is the pressure of the hypothetical cluster defined such as pos- 
sessing the bulk property and fillin g in side the surface of tension. Its deriva- 



tion was given in the literatures |15|, [16|, [17|, lla, ll9| • The present authors have 



given a transparent explanation for the volume term W vo i = —(pp — p a )Vp 
through a grand potential formalism recently [20[]. That is, one can readily 



understand the form of W vo i on the basis of the fact that the reversible work 
of formation of the critical nucleus is the grand potential difference (recall 
that the grand potential 0, is equal to —pV). Nishioka and Kusaka [2l| 
found out that in the case that the /3 phase is incompressible, the volume 
term W vo i = — (pp — Pa)Vp can be rewritten in the form nA/i. That is, they 
integrated (dp,/dp)T = v, which is nothing other than Gibbs-Duhem relation 
for the isothermal case, for the /3 phase (the same procedure was followed 
by Debenedetti and Reiss [181]) . Unfortunately, they concluded incorrectly 
that the form of Eq. (JT|) was valid only for the case of the incompressible /3 
phase such as a nucleation of an incompressible liquid phase in a vapor phase. 
One of the present authors has integrated (dp,/dp)T = v for the a phase in 
case that the a phase is incompressible to get the form of Eq. (pQ) [22]. For 
example, this condition is valid for a bubble nucleation in an incompressible 
liquid phase. For the case of the incompressible /3 phase, we have 

Vr 

W vol = —\j*i){T,p a )-np(T,p0)] 

= —\fip(T,p a )-fi a (T,p a )], (3) 

where T is the temperature, which is assumed to be uniform throughout the 
system. To reach to last expression, we have used the fact that the chemical 
potential of the nucleating phase fjLp(T,pp) is equal to that of the parent 
phase Ha{T,p a ) (that is, the chemical potential is uniform throughout the 
system). For the case of the incompressible a phase, we have 

Vr 

W vo i = — \p a (T,p a ) - (J, a (T,pp)\ 

v a 

Vr 

= — lnp(T,p p ) - n a {T,pp)}. (4) 



v, 



Q 



For extension to the munticoponent system, we may merely follow Nish- 
ioka and Kusaka 21|. In this paper, following Nishioka and Kusaka [21], we 



will expand the chemical potentials for the case that no volume change is 
associated with the a- (3 phase transition. We limit ourselves to the single 
component system for simplicity. 

In Eqs . P) and flU) , fip (T, p a ) - y, a (T, p a ) and fj,p(T,p p )- \x a {T,p p ) are not 
the measurable quantity experimentally directly — at least, the equation of 
state must be measured to be integrated. The "undercooling" (driving force) 



is defined as the difference between targeted chemical potential and that 
at the a-fi equilibrium. However, this terminology is somewhat confusing, 
because a temperature decrement is imagined thereby. Let us use, instead, 
the chemical potential difference in this paper. For a critical nucleus, as 
mentioned above, fj,p(T,pp) and /i Q (T,p a ) is equal with each other and this 
value is the chemical potential of the reservoir. For latter convenience, let 
us define the chemical potential of the reservoir /i re , which is the chemical 
potential of the real system. Also we define the chemical potential at the 
a-j3 phase equilibrium, \i eq . We define the chemical potential difference 

A/i = fi re - fi eq . (5) 

In the case that the a phase is a rarefied gas, for example, this quantity is 
expressed as 

Ap = k B TlnS. (6) 

where S = p a /p eq with p eq being the pressure at the a- (3 phase equilib- 
rium. Some ones inappropriately state that fjLp(T, p a ) — fi a {T,p a ) is mea- 
surable; this is not entirely correct - - this quantity can be expressed as 

/J,p(T,p a ) - H a {T,p a ) = \jMp(T,p a ) - lJteq]-\pl>a(T,P a ) ~ l^eq] = Vp(Pa ~ Peg) ~ 

kBT\ia(p a /p e q) for the case that the a phase is an ideal gas and the f3 phase is 
incompressible (this is, however, not the present concern). [We note that in 
the crystal nucleation form the melt A/i = AS(T m — T) replaces (AS 1 is the 
entropy of melting and T m the melting temperature).] Correspondingly, mis- 
understandin gly /ifl — //„ in Eq. ([T]) is substituted by Eq. ([2D and sometimes 

0, H 0, SB M, m M 



IforoY 

3 A/i 2 l ' 

is used for the barrier height of nucleation, which is obtained by solving 
dW/dR = for R and inserting the solution into W. This situation indicates 
that expressing W in terms of the true driving force, A/i, is strongly desired. 



The other subject of Nishioka and Kusaka [2l| was formulation of the 



reversible work for non-critical clusters. Nishioka and Mori [30| rewrote 
their formula and gave an approximate concise expression. A little later, 



Debenedetti and Reiss [18j drove a similar expression to that in Ref. [21 
This subject is, however, not the present concern. We note that for the 
system including a non-critical cluster, the chemical potential is no longer 
uniform. 



2. Calculations 

In the above way, one has known under what cases the form of Eq. ([1]) 
is valid; those are not the cases which one intuitively imagines. The purpose 
of this paper is to give a form of the volume term of W in this intuitive case 
(the case that the a- (3 phase transition accompanies no volume change). Let 
us start with the p-T phase diagram. The normal case is that the molecular 
volume of the low-temperature phase (/3 phase) is smaller than that of the 
high-temperature phase (a phase). In this case, the phase boundary is a curve 
with positive slope in the p-T phase diagram. The case that the molecular 
volume of the f3 phase is larger than that of the a phase such as the water- 
ice case is an abnormal case. In this case, the phase boundary has negative 
slope. What we will consider is near the point connecting those two resumes. 
For the case that the molecular volumes of two phases are the same over 
a certain region, because the phase boundary becomes a straight vertical 
line in this region in the p-T phase diagram, the pressure no longer induces 
the phase transition. For the case of water-ice, for example, while near an 
atmospheric pressure the water is denser than the ice, at hight pressure the 
ice becomes denser than the water. The phase diagram is as illustrated in 
Fig. [1] (strictly speaking, this is a mere intuition -- the phase boundary in 
p-T phase diagram bends due to transition between the normal ice and a 
high-pressure one as will be discussed in detail in a latter part of Sec. EJ. 
Hereafter, we restrict ourselves on the temperature T = T^ (we sometimes 
omit T^ itself for brevity). We define p eq as the point at which no volume 
change is associated with the a-j3 phase transition. While in the region 
p > p eg , the phase boundary behaves normally, in the region p < p eq it 
behaves abnormally. In the both regions, the /3 phase is more stable than 
the a phase. This means that the pressure-induced phase transition from the 
a phase into the (3 phase takes place. 

jjl(T\p) around p = p eq is schematically drawn in Fig. [2j Except for just 
on the point p = p eq , the a phase is always metastable phase as mentioned 
above. For thermodynamic consistency the a phase is less compressible than 
the (3 phase {n a < Kp). If one imagines a vapor-liquid case as a typical case, 
this observation is, at a glance, surprising because the valor phase is over- 
whelmingly compressible. Considering, for example, a crystalline nucleation 
from a (well compressed) undercooled liquid or a glass phase, this situation 
may be reasonable; focusing on the free volumes of two phases, the crystalline 
solid phase may be more compressible. For the vapor-liquid case, the situa- 



tion illustrated in Fig. [TJ is impossible. Let us express the chemical potential 
in a series expansion around p eq in the pressure up to second order: 



/i(p) - fi 



"i 



' eq 



(P-Peq) ~ ^{P~P 



eqj 



(8) 



where v eq is the molecular volume at the a-(3 phase transition and k the 
isothermal compressibility —(l/v)(dv/dp)T- This equation can be obtained 
by integrating k = —(l/v)(dv/dp)T twice with help of v — (dfi/dp)T as if k 
is constant and then expanding. We obtain from Eq. (JSJ) the series expansion 
solution to p — p eq as 



- eq 



P~Peq = < [M - fi(Peq)] 



K 



2v 



HP) - V{Peq)\ 



rq 



(9) 



We note here that because Eq. (jBJ) is a quadratic equation, we can solve 
it using the quadratic formula and then expand the solution assuming the 
smallness of k to obtain Eq. (Q [one can have a solution x = —c/b — ac 2 /b 3 + 
0(a 2 ) to ax 2 + bx + c = by expanding x = (— b + a/6 2 — 4ac)/2a into 
a power series in a.]. Also, we can get Eq. (jUJ) by solving the quadratic 
equation iteratively assuming the smallness of the coefficient of the second 
order term [one can rewrite the quadratic equation ax 2 + bx + c = into 
x = —c/b — (a/b)x 2 and then the Oth order solution x^ = —c/b and the 1st 
order correction x^ = —(a/b)(x^) 2 = —ac 2 /b 3 ]. Applying Eq. ([9]) with n a 
and Kp begin the isothermal compressibilities at p — p eq for both phases, we 
have 



1 

P/3 -Peq = 

Veq 



Veq) + 



2v 



Or 
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Pa Peq 



■ eq 
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\Vre Veq) ~r 7> \Vre fJ> 



eq i 



2v 



eq i 



Subtracting Eq. ( TTTjl from Eq. 



Pc 



TO]) , we have 

Kp — « c 



2v 2 

eq 



"Or 
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eq J 



(10) 

(11) 

(12) 



Here, we note that if we neglect the compressibility effect, i.e., in the first 
order in //, we have a result that no volume change, no pressure-induced 
phase transition. Putting Eq. ( fl2|) into Eq. (J2J), the reversible work becomes 



W 



Kfi 



Kr 



2v 2 

u eq 



(/Ire ~ VeqfVp + jA. 



(13) 



This equation is, in appearance, of a form far from the commonly used for- 
mula [Eq. CQ)]. 

In this way, we have shown that for the case of no volume change at 
the phase transition, for which the commonly used formula of Eq. ([1]) seems 
intuitively to hold, the reversible work of formation of the critical nucleus 
takes a form far from Eq. ([T]). The volume term of the reversible work of 
formation of the critical nucleus possesses a form including a second order 
in A/i = /i re — [leg, not a first order. The reason why the second order term 
appears is because the second order term is the lowest order term necessary 
for the self-consistent treatment. 

3. Discussions 

Let us discuss a relation of form of Eq. ([3} and the present results 
[Eq. (fl3|) ]. Assuming that the deviation from the equilibrium is small, one 
can make the following expansion. 

PfS ~ Pa = (P/3 - Peg) ~ (p a ~ Peg) 

'' '''' ! ' [ftflM - Wl(Peq)] 



dpp \ 
d/j, J T 



Vp{Pa)=V.eq 



WM ~ t*p(Pa)] + h-O. 



Veq 



\Ma(jPa) ~ ^piPa)] + h.O., (14) 



Veq 



with h.o. standing for higher order terms. Accordingly, instead of Eq. ([3j) 
we have 



W vot = — \p,p(T,p a ) - fj, a (T,p a )] + h.o. (15) 

Veq 

Applying Eq. (jUJ) we can also expand fipiT, p a ) — fi a (T,p a ) as 



v 

fJ,p(T,p a ) - [i a (T,p a ) = -k-fap- K a )(Pa-Peq) 2 - (16) 

Using this equation and Eq. ffTTj) we rewrite Eq. (IT51) into Eq. ffT3l) . In this 
way, we have an additional proof. In addition, versatility of the form of 
Eq. (jHJ) has been shown. 



Here, let us discuss about the reality of the phase transition without 
volume change. In Sec. [2] we have mentioned about the water-ice case. In an 
early phase diagram in p-T space, the phase boundary between the ice III 
and water has been drawn as vertical (see, e.g. Fig. 5.15 of |3l|). In updated 
ones such as in [32|, |33|, |3J, |35|, |36| that boundary is not vertical. Possibilities 
of no volume change, however, remain for ice VII- ice VIII boundary such 
as drawn in the phase diagram 33|, |35|, |36( and for ice Ill-ice IX transition 
regarding bond ordering [37|. There are other possibilities regarding ice- ice 
phase boundaries which are less revealed. A problem about those possibilities 
is that the both phases are solid states so that the phase transitions are slow. 
One may, however, not be disappointed. There exists certainly a point where 
no volume chan ge i s associated with the phase transition on the melting curve 



of graphite [38|, |39|, |40|, |41[. A note is that the slope of the melting curve is 
positive at low pressure and then becomes negative as the pressure increases. 
Therefore, the shape of the melting curve is convex to the right as opposed 
to Fig. [TJ The fundamental framework of the theory is, however, applicable. 



the volume term is of the second order in 
After Nishioka and Kusaka 21 , little of 



4. Concluding remarks 

This paper has made a new finding; a note should be added to some 
textbooks thereby. Intuitively one regards that if the volume change at the 
phase transition can be neglected, the volume term of the reversible work of 
forming the critical nucleus becomes of the form n(/xp — /x a ). This intuition 
should be amended. In this case, 
/x re — Heq, not of the first order, 
concerning description of textbook has been rewritten. The authors hope 
that repeating of the same situation regarding the present finding should 
be avoided. Concretely speaking, the exact form of A/x should be clearly 
written. Without doing so, the misunderstanding that the form of Eq. ([!]) 
would hold for the case of a phase transition without volume change may 
happen. 

The present result may affect the application of the nucleation theorem 

42j, |43|, |44j because the original form is based on the form of nA/x. Although 

the A/x is exactly given by /xp(T,p a ) — ix a (T,p a ) and the volume term is 

given by Eg. (IT51), p ractically d In J/9 In S is calculated such as in literatures 



26l . |45l l46l . 1471 . |48|. Therefore, reconsideration of the application of the 



nucleation theorem is necessary. This is one of the future researches. 



In this paper, we have treated the special case that the phase transition 
accompanies no volume change. Detecting an anomalous behavior in reality 
is one of the future researches. Consideration on the expression in terms of 
l^re — Heq for general cases is in progress. 
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Figure 1: An illustration of the p-T phase diagram near the point where no volume change 
is associated with the a-/3 phase transition. 
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Figure 2: /i-p diagram around p = p eq at T = T>. /i-p relations for both p a > p eq and 
Pa < Peq are schematically illustrated. 
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